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Evidence of accretion triggered oscillations in the 
pre-main-sequence interacting binary AK Sco 
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ABSTRACT 

Pre-main sequence (PMS) binaries are surrounded by circumbinary disks from which 
matter falls onto both components. The material dragged from the circumbinary disk 
flows onto each star through independent streams channelled by the variable gravita- 
tional field. The action of the bar-like potential is most prominent in high eccentricity 
systems made of two equal mass stars. AK Sco is a unique PMS system composed of 
two F5 stars in an orbit with e=0.47. Henceforth, it is an ideal laboratory to study 
matter infall in binaries and its role in orbit circularization. In this letter, we report 
the detection of a 1.3mHz ultra low frequency oscillation in the ultraviolet light curve 
at periastron passage. This oscillation last 7 ks being most likely fed by the gravita- 
tional energy released when the streams tails spiralling onto each star get in contact 
at periastron passage enhancing the accretion flow; this unveils a new mechanism 
for angular momentum loss during pre-main sequence evolution and a new type of 
interacting binary. 

Key words: binaries: close; binaries: spectroscopic; stars: pre-main-sequence; stars: 
individual: AK Sco ; circumstellar matter, accretion: accretion disks 
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1 INTRODUCTION 

The formation of binary stellar systems is not yet fully un- 
derstood (Pringle 1991). A key problem is the abundance of 
close binaries that implies the action of efficient mechanisms 
for gravitational energy loss and angular momentum trans- 
port further than the tidal interaction between the stars. 
Star-disk interaction during the PMS evolution is thought 
to play a crucial role (Clarke & Pringle, 1991). Accretion 
disks can either take up or release angular momentum in 
the system that works like a multi-body gravitational sys- 
tem (Artymowicz & Lubow, 1994; Bate & Bonnell, 1997; 
Hanawa et al 2010). The details depend on the mass ratio 
between the two stars and on the orbit eccentricity. Highly 
eccentric orbits favour the formation of spiral waves within 
the inner disk that do channel the flow as the accreting gas 
streams onto each star. 

AK Sco stands out among the PMS binaries by the 
high eccentricity of its orbit, e = 0.47; the two F5-type 
stars get as close as 11.3 stellar radii at periastron pas- 
sage (see Table 1, for AK Sco main parameters). UV ob- 
servations carried out with the Hubble Space Telescope 
(HST) unveiled the presence of a dense and extended struc- 
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ture around the stars with thermal properties alike those 
of the atmospheres/magnetospheres of late type PMS stars 
(Gomez de Castro 2009). The structure extends up to 1.7 R, 
around each component and it is disturbed by a velocity field 
(a ~ 100 km s" 1 ) that exceeds significantly the sound speed 
in the UV radiating plasma, 30.5 km s - . Thus, HST obser- 
vations indicate that there are unresolved macroscopic mo- 
tions in the dense and warm circumstellar material. Given 
the eccentricity of the orbit, we expected mass transfer to 
be significantly enhanced at periastron passage, forcing the 
response of the stellar magnetosphere to the accretion flow. 

During PMS evolution, stellar magnetospheres play a 
key role as dissipative interfaces between the star and the 
disk. They do absorb and reprocess part of the angular mo- 
mentum excess of the infalling material and channel it into 
open holes in the stellar magnetic configuration. Henceforth, 
magnetospheric radiation (ultraviolet and X-ray) is a very 
sensitive tracer of accretion. For this reason, we monitored 
AK Sco with the XMM-Newton space telescope at three 
different phases, including the periastron passage (more de- 
tails on the monitoring campaign can be found in Gomez de 
Castro et al. 2012). However, it was only at periastron pas- 
sage that an oscillation was detected at the rise of the light 
curve. The observations and data analysis are described in 
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Table 1. Main parameters of AK Sco 



Property 


Value 


Reference 


Projected semi-major axis 


asini = 30.77±0.12 R Q 


Andersen et al 1989 


Eccentricity 


e= 0.47 


Andersen et al 1989, Alencar et al. 2003 


Orbital period 


P=13.609 d 


Andersen et al 1989, Alencar et al. 2003 


Inclination 


i = 65° - 70° 


Alencar et al. 2003 


Age 


10-30 Myrs 


Alencar et al. 2003 


Spectral type 


F5 


Alencar et al. 2003 


Stellar Mass 


M, = 1.35 ± O.O7M 


Alencar et al. 2003 


Radius 


R* = (1.59 ± 0.35).% 


Alencar et al. 2003 


Projected rotation velocity 


v sini = 18.5 ± 1.0 km s~ 


1 Alencar ct al. 2003 


Bolomctric flux 


6.33 X 10~ 9 erg s^cm -2 


Andersen et al 1989 


Extinction: Ay 


0.5 mag 


Manset et al. 2005 


Extinction: R 


4.3 


Manset et al. 2005 


Parallax 


9.73 mas 


Van Leeuwen 2007 



Section 2. In Section 3, the possible sources of the oscillation 
are analysed. A brief summary is provided in Section 4. 



2 OBSERVATIONS AND DATA ANALYSIS 

AK Sco observations were obtained on March 15th, 18th 
and 22nd, 2011. They were performed at phases 0.0, 0.15 
and 0.48, corresponding to observation identifications, ID 
0651870201, ID 0651870301 and ID 0651870401, respec- 
tively. Phase 0.0 corresponds to periastron passage. The ex- 
posure time of each observation was approximately 25 ks. In 
each observation, both the OM and the EPIC instruments 
were used to study the ultraviolet (UV) and X-ray evolution 
of the magnetospheric plasma, respectively. 

The OM was used in Image and Fast modes with the 
UVM2 filter (2310 A). The AK Sco spectrum is dominated 
by magnetospheric radiation (mainly the Balmer contin- 
uum and the Fe II resonance UV multiplets) in the wave- 
length range of this filter (2000 A - 2700 A) thus the photo- 
spheric contribution is negligible (Gomez de Castro & Fran- 
queira 1997). Each observation consisted of four exposures of 
4400 s. Data were processed with the XMM-Newton Science 
Analysis System. We obtained light curves with a time reso- 
lution of one second, re-binned to 100 s to increase the signal 
to noise ratio. Fast oscillations were detected by OM only 
from 4000 s to 12000 s during the observation performed 
at periastron passage; notice the three clear peaks between 
6 ks and 9 ks in Figure 1. The peaks have the same struc- 
ture with a first heading maximum and a second weaker 
peak. The strength of the main peaks is above 2o in the 
light curve. The X-ray signal was too low to detect similar 
oscillations in that band. 

To analyse in more detail the oscillations, we used the 
method applied by Torrence & Compo (1998) to the "El 
Nino" oscillation in the Atlantic Ocean. The analysis is 
based in the wavelet transform technique, a modified version 
of the Fourier analysis for the detection of quasi periodic sig- 
nals. This method is essentially a form of windowed Fourier 
transform; it consists of convolving the original signal with 
a wavelet (or mother) function that is localised in both time 
and frequency domain. Examples of the use of this method 
in Astrophysics can be found in Bijaoui 2004, Welsh et al. 
2006 or Mathioudakis et al 2006. 
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Figure 1. Upper panel: OM light curve of AK Sco at periastron. 
Time binning is 100s . Lower panel: wavelet power spectrum of 
the upper light curve. The Y-axis is the period (P = , with 

v being the frequency of the power spectrum) . Contour levels are 
represented for la and 2a significance levels. The dashed area 
represents periods for which results are not statistically reliable. 



We used the Morlet function as mother function (Mitra- 
Kraev et al. 2005), since it is especially well suited to detect 
non- stationary oscillations (Torrence & Compo, 1998). The 
resultant power spectrum is also shown in Figure 1. The os- 
cillation period was found to be 790^0 s, independently of 
the temporal sampling of the OM light curve. The amplitude 
is 1.5±0.2 counts/s; about a 6% oscillation over a 25 count/s 
signal level. In addition, we reconstructed the light curve by 
frequency domains. The frequency-filtered light curves are 
shown in Figure 2 and the sum of the three curves gives the 
observed OM light curve in Fig. 1. Retaining only the low 
frequency (large periods) permits to observe the large scale 
shape of the detected signal. Notice that the oscillation oc- 
curs in the middle of the rise. The middle panel shows the 
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Figure 2. Reconstructed light curve after frequency filtering. The 
periods for the filtering are written in each panel. The scales are 
normalised for convenience. The low frequency (large periods) is 
shown in the top panel. The middle panel shows the reconstructed 
light curve for frequencies around the peak of the wavelet power 
spectrum (see Fig.l); the dashed lines mark the +/-2<r rms of the 
light curve. The reconstructed light curve for high frequencies is 
shown in the bottom panel. 

reconstructed light curve for frequencies around the peak 
of the wavelet power spectrum (see Fig.l). This light curve 
shows a clear modulation. For high frequencies, the recon- 
structed light curve shows basically white noise. 



3 THE SOURCE OF THE OSCILLATION 

This ultra low frequency (ULF) oscillation represents a long 
wavelength mode excited around periastron passage and 
lasting only for few oscillations at the rise of the UV light 
curve. 

At periastron, the distance between both stars shrinks 
to 11.3 R* from a maximum of 31.4 R* distance at apas- 
tron passage. Studies of the structure and evolution of cir- 
cumbinary disks predict the formation of circumstellar disks 
around each component (see e.g. Hanawa et al 2010). Cir- 
cumstellar disks are constrained within the Roche lobe thus 
matter falling from the circumbinary disk is being pulling 
towards the Lagrangian points LI, L2 and L3 feeding them. 
Matter orbits in the same direction in both circumstellar 
disks thus, when disks get in contact at periastron passage, 
a shear instability could be produced at the interface driving 
to angular momentum annihilation. An accretion outburst 
should be produced at this time. As shown in the light curve 
in Fig.s 1 and 2, the UV flux increases just before the peri- 
astron passage and the oscillation is observed in the rising 
of the light curve. 

However, the precise mechanism producing the 790 s 
oscillation is unclear. One might expect that the ULF os- 
cillation is tracing thermal oscillations in the inner border 



of the circumstellar disks. The avalanche of gas from the 
contact area between both disks, i.e. from the outer radius, 
would increase the density and opacity at the inner border 
of the circumstellar disks driving to thermal oscillations, as 
often observed prior to the outburst in interacting binaries 
(see e.g. Gaesincke et al. 2006). However, the ULF oscilla- 
tion period is significantly larger than the thermal timescale, 
T t h ~ 1/tvfl with a the parameter prescribing transport in 
a thin circumstellar accretion disk and ft the Keplerian fre- 
quency at the inner border of the disk. 

This period, however, is close to that observed in the 
magnetosonic oscillations of flux tubes in the Sun (de Moor- 
tel 2009). Periods range from 284 s in coronal loops (McE- 
wan & de Moortel, 2006) to 600-900 s in coronal plumes (Of- 
man et al. 1999). They are thought to be triggered by oscil- 
lations in the loop foot-points that can either be impulsive 
(Nakariakov et al. 2004) or associated with solar p-modes 
(McEwan & de Moortel, 2006). Magnetoacustic oscillations 
with similar periods have only been detected in the very 
magnetically active star AT Mic apart from the Sun (Mitra- 
Kraev et al. 2005). But AK Sco cannot be considered as a 
very active star; it is an F5 star with very moderate X-ray 
flux. Henceforth, only a powerful external trigger could ex- 
cite these oscillations to the level of making them detectable 
over its prominent UV excess. Only, the rapid enhancement 
of the accretion rate at periastron could play that role. 

During PMS evolution stellar magnetospheres extend 
to the inner border of the accretion disk since they rotate 
more slowly than the Keplerian disk. Infalling matter piles 
up at the magnetosphere boundary while its angular mo- 
mentum excess is transferred to other scales or dissipated 
(Ghosh &: Lamb 1979, Goodson et al 1999, Romanova et 
al 2012). Magnetic flux tubes are expected to connect the 
stellar surface with this sheared interface layer (Favata et 
al. 2005). Rapid plasma penetration at the base of the loop, 
in the inner border of the accretion disk, may excite oblique 
magnetosonic shock waves producing an impulsive flare. The 
additional energy filling the magnetic flux tube with hot 
plasma determines the relative amplitude of the oscillation 
thus AI/I — (cs/va) 2 = 0.06, with c s the sound speed 
and va the Alfven speed in the loop. Assuming that the 
thermal conditions within the loop are similar to those de- 
termined for the average corona from the X-ray spectrum, 
the magnetic field in the loop is estimated to be B = 52.8 G, 
a modest value as expected for an F5 type star. However, 
the natural oscillation period of the magnetosonic modes 
is: r = 2L/jcr, with L the length of the loop, j the order 
of the harmonic and (1/ct) 2 = (1/cs) 2 + (1/va) 2 - Under 
coronal conditions, r = L/cs, for the first order mode. This 
corresponds to a natural loop length of 0.14R* for AK Sco, 
that seems to be very small for loops connecting the star to 
the inner border of its circumstellar disk. Non to mention, 
the high degree of synchronisation required to observe the 
oscillation simultaneously in both stars. 

The spiky nature of the oscillations suggest that they 
could be associated with bursts, typically triggered by recon- 
nection events. However, the large amplitude of the outburst 
(6%) indicates that there is an external trigger as accretion. 
In PMS stars, the stellar field is assumed to be anchored 
in the inner part of the disk creating a sheared layer be- 
tween the rigid body rotation of the star and the Keple- 
rian rotation of the disk. This interaction has a profound 
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Figure 3. In PMS stars, the interaction between the stellar mag- 
netic field and the disk twists the stellar field lines due to the 
differential rotation. The toroidal magnetic field generated out of 
the poloidal flux and the associated pressure tend to push the 
field lines outwards, inflating them, and eventually braking the 
magnetic link between the star and the disk (boundary between 
regions I and II). Three basic regions can be defined: Region I 
dominated by the stellar wind, Region II dominated by the disk 
wind and Region III dominated by stellar magnctosphcric phe- 
nomena. The dashed line traces the boundaries between this three 
regions. The continuous lines indicate the topology of the field and 
the shadowed areas represent regions where magnetic rcconncc- 
tion events are likely to occur, producing high energy radiation 
and particles (from Gomez dc Castro 2004). 

influence on the star and the accretion flow but also acts 
as a dynamo that transform part of the angular momen- 
tum excess in the inner disk into magnetic field amplifica- 
tion that self-regulates through quiescent periods of field 
building-up and eruptions when the energy excess is released 
by magnetic reconnection in the current layer (see Figure 
3). There are however, a great number of uncertainties in 
the way the system self-regulates that, in turn, depends on 
magnetic diffusivity (see for instance, Goodson & Winglee, 
1999). Typical diffussivities range from r\ — 92.34 m 2 s _1 in 
the dense magnetospheric components radiating the ultra- 
violet to r\ = 0.09 m 2 s _1 in the low density, X-ray radiating 
plasma. If the ULF oscillation was associated with quasi- 
periodic oscillations in reconnection events within the disk- 
star boundary layer, the inferred physical scales, L , would 
be small since Lo = (77790s) 1 / 2 and compatible with the 
expected thickness of the layer. 



4 CONCLUSIONS 

We report the detection of an ULF oscillation in the UV 
light curve of the PMS binary AK Sco at the rise of the light 
curve. This is the first detection of an accretion triggered os- 
cillation in a PMS binary and makes AK Sco an ideal labora- 
tory to study angular momentum evolution in these systems 
as well as close binary formation. This work unveils a new 



type of interactive binary where outburst are driven by an- 
gular momentum loss between circumstellar disks. However, 
further observations are required to identify the physical na- 
ture of the oscillation. The frequency, timing and strength of 
the oscillation at various wavelengths and thermal regimes, 
needs to be investigated, as well as the repeatability of the 
phenomenon. The detection of the ULF oscillation opens up 
the way for detailed studies on the properties of PMS mag- 
netospheres and their interaction with the accretion flow. 
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